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To investigate the large scale structure of the universe from a field-theoretical perspective, and to give a 
unifying understanding of several seemingly-unrelated observed features of clusterings of galactic objects, such 
as the dependence of clustering amplitude on mass and luminosity, the departure of the correlation function ^(r) 
from the power-law oc (rg/ry '' , the oscillatory behavior of ^cc(r) for clusters on very large scales, and the scaling 
of ro with the mean separation d of clusters. 

We present a coarse-grained field theory of density fluctuations for a Newtonian self-gravitating many-body 
system, and apply it to a homogeneous Universe with small density fluctuations. The theory uniformly treats 
the clustering of all types of galactic objects in terms of the field of density fluctuations. 
I The Jeans length Aq, a unique physical scale for a gravitating system, appears naturally as the characteristic 

^— ■») , scale underlying the large scale structure. Under Gaussian approximation the analytic expressions of ^(r) and 

' P(k) are obtained. The correlation amplitude is proportional to the galactic mass, ^(r) oc m, and ^(r) departures 

' from the power-law (ro/r)''^ around r ~ 20 Mpc, and is oscillating over large scales ~ 100 /i^'Mpc and 

^ , damped to zero . The spectrum amplitude is inversely proportional to the galactic number density, P(k) oc 

^ ■ and P(k) oc k^--^ in the range k ^ (0.05 ~ 0.3) hMpc"'. The scaling ro oc ^/O-^-O-^ holds quantitatively. 

' These preliminary analytic results under Gaussian approximation already qualitatively explain these pro- 

nounced features of the observed large scale structures. Yet on small scales the predicted clustering is insuffi- 
cient, and the peak of P(k) is too sharp. These deficiencies are due to neglecting higher order nonlinear effects, 
which need to study further. 
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O ■ I. INTRODUCTION 

^! In past decades, great progresses have been made in understanding the the large scale structure of the 
O universe. Not only observations of the major galaxy surveys such as SDSS, 2dF, APM, et. al., (Tegmark 
^ 2004; Zehavi 2005; CoUess 2001; Hawkins 2003; Madgwick 2003; Percival 2005; Maddox 1996; Padilla 
^ 2003) have revealed the cosmic structures of increasingly large sky dimension, theoretical studies have also 
achieved important results, through both analytical modelling and numerical simulations. The correlation 
• ^ functions of galactic objects, are key physical quantities, which contain important information, both dy- 
^ namic and statistical, of the cosmic structure. Observations have shown that the correlation functions of 
various galactic objects, such as galaxies, groups of galaxies, quasars, and clusters, are usually fitted by a 
single canonical power-law ^(r) = (ro/ry in the range r < 4ro with the index y ~ 1.7 and the 'correlation 
scale' ro ~ 5/z~'Mpc for galaxies, ~ 10/?"^Mpc for galaxy groups, ~ 20/7"'Mpc for clusters, respectively 
(Peebles 1993; Bahcall 1996; Bahcall et al 2003). Also ro increases with mass, richness, and luminosity 
of galactic objects. Moreover, ro scales with the mean cluster separation d as ro{d) oc (Bahcall 
1996; Croft 1997; Gonzalez 2002; Bahcall et al 2003; Zandivarez, Merchan, Padilla 2003). On large scales 
^(r) departures from the power-law, drops steeply to zero, and seems to have an oscillating behavior on 
very large scales > 100/j^'Mpc. These observed features have been thought to be a theoretical challenge 
(Bahcall 1996; Bahcall et al 2003). From view point of dynamics, the Universe filled with various galactic 
objects is a many -body system with attractive gravity between bodies. A powerful, systematic approach 
to deal with statistical dynamics of this kind of systems is through introducing the coarse-graining order 
parameter, which converts the degrees of freedom of discrete particles into the effective continuous field. 
Therefore, a field theory arises, which combines both the statistical and dynamical properties of the system. 
In particular, the method has been very successful in revealing the correlation properties of fluctuations 
for a variety of homogeneous systems near phase transition (Goldenfeld 1992). Therefore, it should be 
appropriate to employ the method to study the density fluctuations of the homogeneous universe as a self- 
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gravitating systems. In this paper we will formulate such a field theory of density fluctuations and apply it 
to the Universe on large scales. It turns out that our analytic result at the level of Gaussian approximation 
already interprets qualitatively the above-mentioned several features simultaneously. Different from the 
conventional analytical studies, such as the perturbation theory (Bernardeau 2001) and the thermodynamic 
formulation (Saslaw 2000), our approach is a field-theoretical study, like those used in condensed matter 
and in particle physics. For instance, in our theory, the 2-point correlation function turns out to be the 
Green's function, and the differential equation that it satisfies can derived explicitly. 



II. FIELD EQUATION OF CORRELATION FUNCTION 



Consider a system of galaxies of mass m in the Universe with the Hamiltonian H = 2, p J2m 



Yu^j Girp- /\ri - Tjl. The grand partition function of this system at temperature T is 



where z is the fugacity. By using the Hubbard-Stratonovich transformation (Zinn- Justin 1996) it can be 
cast into a path integral: Z = J Dcpexp [-/3~^ j d^r £.((/))], where is the gravitational field, /3 = AnGmjc] 
of dimension m~\ = {T ImY^^ the sound speed, the effective Lagrangian is X(0) = ^(V0)^ - fc^e'^, 
and kj = (AnGpo/c^^^^ is the Jeans wavenumber. The term -k^e'^ has a minus sign because gravity is 
attractive. By the Poisson's equation v^0(r) -i- kje'^ = the mass density p is related to the field by 
p(r) = mn{r) = poe'^'-^K So po is the constant mass density when = 0. We define a dimensionless re-scaled 
mass density field 

0(r)^/«=p(r)/po, (2) 

viewed as the coarse-grained order parameter, and introduce a Schwinger type of external source J coupling 
with ij/ in the effective Lagrangian 

X(<A,/) = ^(yj -k]iff-Jif,. (3) 

The source 7 is used to handle the functional derivatives with ease. The equation for i/r in the presence of 7 
is 

v2 - i(V(/^)2 + kW + Jif? = 0, (4) 



the generating functional is 



Z[J] = J Di/,cxp[-/3'' Jdh£{i/f,J)l (5) 



the n-point connected correlation function is 

G("^(ri,...,r„) = <5(A(ri)...50(r„)) 

„^MnZ[7^ (6) 
^ 57(ri)...57(r„)" 

where 6iff(r) = ij/ir) - (i/^(r)) is the fluctuation field about the expectation value (<A(r)). The standard way to 
obtain the equation for the 2-point correlation function Gf\v) is the following (Goldenfeld 1992). Taking 
expectation value of Eq.® in the presence of 7, using the approximation {(yij/Y - (v(i/r)/)^/(i/r)y. 
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which is equivalent to neglecting high order terms of 5ip, taking functional derivative of this equation w.r.t. 
the source /, and setting 7 = 0, we end up with 

Gf\r)+ (^)2Gf (r) - A v • siGf{r) 

w) w) 

+2k]{ilf)Gf{r) = -<iA)W^(r), (7) 

where Gf\r - r') = /3d{i^(r)) j=o/dJ(r') has been used, and (i/r) = (i//) is the background field depending 
on the position r in general. We remark that when higher order terms in 1 /i/^ =^ 1 /i/tq + ^i/^/i/^q + ... were 
kept, one would end up with nonlinear equation much more complected than Eq.©. Assume now the 
homogeneity of the Universe on large scales. Then the background density field is constant {if/) = if/o = I. 
Then Eq.© reduces to the Helmholtz equation with a point source: 



AnGm 



^(r) + Ik'j^ = —S'^'ir), (8) 

where we have used a simple notation ^(r) = Gp^(r). The approximation level of Eq.® is essentially 
that of the Gaussian model in the Landau-Ginzburg theory(Goldenfeld 1992). The solution is the Green's 
function 

/ cos(fcor) sin(fcor) \ 
^(r)=A\c + (l-c) — I, (9) 

where the amplitude A = Gm/cl, the parameter c determines the contributions from independent cos and 

sin modes, and the characteristic wavenumber = ^kj. Thus the wavelength Aq = Aj/ V2 occurs as the 
characteristic scale that underlies the large scale structure, where Aj is the Jeans length. The 'correlation 
scale' ro in the canonical power-law is only empiric not fundamental. The amplitude A oc Gmlc]rQ, the 
ratio of the gravitational correlation energy Grn^/ro to twice the kinetic energy mc^ of peculiar motion of 
galaxies. So a greater gravitational energy yields a larger correlation amplitude, similar to the conclusion 
by Saslaw (Saslaw 2000). A comparison with the Landau-Ginzburg theory will be suggestive. Recall that 
for a fluid, with the dominant force between molecules being repulsive, the Gaussian approximation yields 
a correlation function oc e~''^^/r, where rj is the correlation length of phase transition (Goldenfeld 1992). 
Contrarily, Eq.® is oscillating, instead of decaying, because gravity is attractive. In this sense, the self- 
gravitating system is always at the critical point of 'phase transition' with 77 ^ oo (Saslaw 2000). Therefore, 
it is pertinent to study the density fluctuations of the homogeneous universe within the framework of the 
field theory, where the coarse-grained order parameter is the rescaled density field if/. If the gravity were 
repulsive, the term Ikjifr in Eq.® would have a minus sign, and Eq.® would be of the same form ^(r) oc 
Q-Mr fQj. ^ fluid. 

^{r){r - 

r^T^'^rdr, and Eq.®, one gets Wp(rp) = 7rA[-cA^o(^o''p) + (1 - c)Jo{kQrp)], where A^o is the Neumann 
function and Jq the Bessel function. So the amplitude is also proportional to the mass, Wp oc m, and 
Wp(rp) has a strong departure from the power-law and drops steeply to zero, agreeing with the observations 
(Budavari 2003, Zehavi et al 2002; 2005). 



III. IMPLICATIONS 



Now examine some implications of ^(r) for the large scale structure. Although, as the Gaussian ap- 
proximation, Eq.® does not yet give a detailed account of of large scale structure, they already reveal the 
following pronounced features. 

1) Galactic objects in our theory are all treated as gravitating particles, and have the same form of cor- 
relation functions as Eq.®, diff"ering only in A and ^o- Indeed, the observed ^(r) for galaxies, groups, and 
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clusters all have a similar shape, for instance, the correlation function of RELEX clusters divided by ~ 10 
coincide with that of EPS galaxies (Guzzo 2000). 

2) Since the amplitude A oc m, galactic objects with a greater mass have a higher correlation. Be- 
sides, since richness is proportional to m by definition and more luminous galaxies generally have greater 
mass, this result naturally explains why the correlation amplitudes increase with richness and luminosity of 
galactic objects without recourse to the bias model. Attempts have been made to explain this in terms of the 
statistics of rare peak events (Kaiser 1984). Moreover, early type galaxies (ellipticals) have a mean num- 
ber density n lower than that of late type (spirals) and are more likely located in denser regions, thus they 
represent a greater mass of the region (Davis, Geller 1976; Madgwick 2003) and have a higher correlation. 

3) ^(r) has a departure from the canonical power law, as is observed. Plot ^(r) in FIGlHto fit with ^gg{r) 
of the 2df galaxies (Madgwick 2003). The parameters A ^ 3 h^^Mpc, c - 1.7 , and fco - 0.05 hMpc"^ 
are taken. It is seen that the departure from the power-law occurs at ~ 8 h 'Mpc, and ^(r) goes to zero at 
- (20 ~ 30) h 'Mpc. But on small scales < 2 h^^Mpc, ^(r) oc r \ lower than the observed. This small-scale 
deficiency of correlation is expected since significant clustering effects of higher order terms of 6>p have 
been neglected in Eq.®. We expect these nonlinear effects will improve this small-scale correlation. For 
groups and clusters, fitting can be done similarly, requiring a larger A and a slightly smaller ^o- However, 
the redshift-space correlation function from the SDSS luminous red galaxies crosses to zero only after 
100/z"^Mpc (Eisenstein et al 2005; Huetsi 2006). It is difficult for our theory, and it might indicate that 
either the universe is still less relaxed, or the higher order nonlinear effects need to be taken into account. 

4) The shape of observed correlation function of galaxies depends on their 'color'. The red galaxies are 
more luminous with a slightly higher mean number density n, and exhibit a steeper and higher-amplitude 
correlation function than the blue galaxies (Zehavi 2002; 2005). This is explained by Eq.®, which, for a 
larger ko oc (nmy^^, is steeper and has a higher amplitude on small scales. 

5) As mentioned at the beginning, the observations show that the 'correlation scale' ro increases with the 
mean separation d in such a way as ro{d) oc (Bahcall 1996; Bahcall et al 2003; Croft 1997). The 
mechanism underlying was thought to be a scaling behavior (Szalay, Schramm 1985). 

In our theory, if we adopt the power-law fitting r^id) oc ^'^'^ with ^ being given in (|9l), then, for typ- 
ical clusters, rQ{d) ~ 20cos^^^-^ [(30/(i)^^^ ], increasing in the range J = 10 ~ 80h"^Mpc, as shown in 
FIGlll Although it is only a approximate result, the trend of the profile is in qualitative agreement with the 
observation. 

6) On very large scales, ^(r) is a damped oscillation about zero. Some surveys have indicated this feature 
with Ao ^ 120 ~ 200 h'^Mpc (Broadhurst et al 1990; Tucker 1997; Einasto et al 1997; 1997b; 2002; 
2002b; Tago 2002; Eisenstein et al 2005; Yanata 2005). Plotting ^(r) in FIGHto fit to the observed ^^dr) 
of clusters requires A =^ 50 h^^Mpc and ko ^ 0.05 h Mpc*^ Thus the Gaussian approximation in Eq.® 
gives a good description of the oscillatory correlation on very large scales. 



IV. POWER SPECTRUM 



The power spectrum is the Fourier transform of the correlation function, and can be also derived directly 
by Fourier transforming Eq.®, 

P(k) = ^-T^- (10) 
2n(f )2 - 1 

Firstly, this predicts that P(k) has a higher amplitude for galactic objects with a smaller spatial number 
density n, agreeing with a whole variety of surveys. For instance, the observed P(k) for all Abell clusters is 
lower than that for Abell clusters with measured redshifts, and the same fact is also found for APM clusters 
(Einasto 2002b). This is easily understood by Eq.lfTOl) as all clusters have a higher n than those with 
measured redshifts. Secondly, P(k) has a sharp peak at k = kg, corresponding to the periodic wavelength 
Aq = 2n/ko. Taking the observed value /lo - 120 ~ 200 h^^Mpc yields the peak location ^ 0.031 ~ 0.052 
hMpc ^ The peak has been, to some degree, indicated by several cluster surveys (Ratcliffe et al 1996; 
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Tadros, Efstathiou, Dalton, 1998; Hoyle et al 1999; Collins 2000; Schuecker et al 2001; Carretti et al 2002; 
Einasto et al 2002; 2002b; Tago et al 2002; Padilla, Baugh 2003; Tegmark 2004; Percival 2005). 

Thirdly, as shown in FIGlH for large scales k =^ (0.05 ~ 0.3) hMpc"\ (flOl) gives P(k) oc k~^-^, agreeing 
well with the major surveys (Peacock 1997; 1999; Carretti et al 2002; Dodelson et al 2002), whereas the 
power-law with y ~ 1.7 would yield oc k~^-^. But for small scales k ^ (0.4 ~ 3) hMpc~\ Eq. (fTOl) gives 
P{k) oc k'^, the observed is oc k~^-^ (Peacock 1999; Dodelson et al 2003). This small-scale deficiency is 
the same phenomenon occurring in ^(r) noticed before, also due to neglect of the nonlinear terms. But 
there are two more difficulties for P(k) in Eq.(fTOl). It takes on negative values on large scales with k < ko, 
and its peak is too sharp. We expect that, when the nonlinear and frictional effects are properly taken into 
account, the negative power at ^ < ko may also be alleviated and the peak may be smoothed to certain 
extent. These are obviously different from the observed power spectrum. A real spectrum can not have 
an infinitely sharp peak, which should be smoothed by some damping effects. Recent observations have 
indicated that the measured P{k) has a rather flat peak with large systematic errors on large scales (Einasto 
et al 2002b; Outram et al 2003; PadiUa, Baugh 2003; Tegmark et al 2004; Mailer et al 2005; Percival 2005). 

V. CONCLUSION 

The effective field theory for density fluctuations presented in this paper has the Jeans scale - (0.04 ~ 
0.06) hMpc ' as the unique scale underlying the large scale structure. Even at the level of Gaussian approx- 
imation its preliminary results already interpret simultaneously several important features of the observed 
large scale structures, some of them have been thought to be puzzling and seemingly-unrelated. But the 
Gaussian approximation does not yet give a sufficient clustering on small scales. Therefore, it is worth 
while to explore further the field theory beyond the Gaussian approximation, and to take into account of 
higher order terms of 6tff neglected in Eq.®. This is expected to enhance the low amplitude of clustering 
on small scales and smooth out the sharpness of peaks in the power spectrum. The study of these nonlinear 
effects is still underway. 
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FIG. 3: The large scale oscillations of ^(r) compared with the observed of clusters (Einasto et al 2002; 2002b). 
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FIG. 4: The analytic P{k) in Eq.l[TOll compared with the SDSS galaxies (Dodelson et al 2002) in the range ^ (0.1 ~ 5)hMpc"'. 
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